QCD traveling waves at non-asymptotic energies 
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Using consistent truncations of the BFKL kernel, we derive analytical traveling-wave solutions 
of the Balitsky-Kovchegov saturation equation for both fixed and running coupling. A universal 
parametrization of the "interior" of the wave front is obtained and compares well with numerical 
simulations of the original Balitsky-Kovchegov equation, even at non- asymptotic energies. Using 
this universal parametrization, we find evidence for a traveling-wave pattern of the dipole amplitude 
determined from the gluon distribution extracted from deep inelastic scattering data. 



I. INTRODUCTION 

In the past decade, there has been a large amount of work devoted to the description and understanding of 
quantum chromodynamics (QCD) in the high-energy /density limit corresponding to saturation. On the theoretical 
side, progress is made in obtaining non-linear QCD equations describing the evolution of scattering amplitudes towards 
saturation. On the phenomenological side, the discovery of geometric scaling Q in deep inelastic scattering (DIS) at 
HERA was a stimulating observation. Indeed geometric scaling has a natural explanation in terms of traveling- wave 
solutions of the Balitsky-Kovchegov (BK) non-linear equation ^i^. 

The BK equation is an equation for the evolution in rapidity Y of the scattering amplitude N{r, b, Y) of a dipole 
(a colorless qq pair) off a given target, r is the transverse size of the dipole and b its impact parameter. The rapidity 
Y is the logarithm of the squared total center-of-mass energy of the collision. Assuming that the amplitude N is 
b— independent, the BK equation has been shown to lie in the same university class as the Fisher-Kolmogorov- 
Petrovsky-Piscounov (FKPP) equation j^, extensively studied in statistical physics. 

The FKPP equation is an evolution equation for the function u{x, t) where x is a space variable and t is the time. It 
admits asymptotic solutions which are functions of the sole variable x—vt, namely traveling-wave solutions u{x—vt) 
where v is the speed of the wave. It comes from a critical regime obtained by the competition of the exponential 
growth in the linear regime and the non-linear damping. This mechanism for the formation of traveling waves is more 
general and applies to the BK equation for which time is replaced by rapidity and the logarithm of the dipole 
size plays the role of the position. It has even been extended beyond the b— independent case to processes with 
non-zero momentum transfer. 

However there are limitations for applying these solutions precisely because they are asymptotic, i.e. they require 
very large values of Y to appear; also analytical expressions are only known for relatively small values of r — |r|. 
This leads to a limitation of analytical predictions to the tail of the wave {N ^1) which is essentially determined 
by the linear regime. For instance, the "interior" of the wave is not reproduced; it corresponds to the transition to 
saturation, a intermediate regime between the tail (7V<C 1) and the saturated region {N ^ 1). That moderates the 
phenomenological impact of those mathematically powerful results. 

In order to circumvent these difficulties, it is necessary to extend the study to the non-asymptotic regime in a 
region of larger N, which implies to use more information on the non-linear term of the equation than the previous 
approaches. For this sake, a new approach to study the BK equation has been developed in In short, the guiding 
idea is to try and find exact solutions of approximate QCD equations rather than approximate solutions of exact QCD 
equations. 

The goal of this paper is to derive these analytical solutions and to develop their applications in various cases 
including the important running-coupling case. This allows to make a direct and model-independent comparison for 
saturation between QCD predictions and phenomenology. 

The plan of the paper is as follows. In Section II, we define consistently truncated BK equations for which we derive 
exact traveling-wave solutions in the fixed-coupling case (HA) and study the validity of the method by comparison 
with numerical simulations of the original BK equation (HB). In Section HI, we find solutions to the BK equation 



t on leave from the fundamental theoretical physics group of the University of Liege. 
^ URA 2306, unite de recherche associee au CNRS. 
*Electronic address: marquetOspht.saclay.cea.fr^ 
^Electronic address: pesch@spht.saclay.cea.fr 
5 Electronic address: gsoyez@spht.saclay.cea.fr 



with running couphng. In Section IV, we apply the results to DIS phenomenology. Conclusions and outlook are given 
in Section V. 



II. THE FIXED-COUPLING CASE 



Let us consider the Fourier transform of the dipole scattering amplitude defined as 

r°° dr 

M{k,Y)^ / ^Jo(fcr)7V(r,r) . (1) 
Jo f 

The BK equation for N{k, Y) then reads 

dYAf = x{-dL)M-N^ (2) 

where L = log [kF' /k^^ and the rapidity Y is measured here in units of q;= ""J^" where the strong coupling constant 
as is kept fixed. 

x(7) = 2V'(l)-^(7)-V'(l-7) (3) 

is the Balitsky-Fadin-Kuraev-Lipatov (BFKL) kernel JJi]. /cq is an arbitrary scale. 

It has been proved that the equation (|2Jl admits asymptotic traveling-wave solutions N{L—v Y) where v is the 
speed of the wave. Starting with the initial condition M{L,Yf)) exp(— 7iL), traveling-wave solutions are formed 
during the y-evolution with the following value for the speed: 

• "pushed front": < 7i < 7c; v = x{li)hi, 

• "pulled front": 7^ > 7^; v ^Vc = x{lc)hc, 

where 7c is solution of the equation x{l)ll — x'il) and is called the critical anomalous dimension. Its value is 
7c = 0.6275 leading to the critical velocity Wc = 4.883. Note that Q only the large— L behavior of the initial condition 
matters to determine whether one lies in a pushed or pulled-front case. 



A. Exact traveling-wave solutions of truncated BK equations 

In order to avoid mathematical difficulties related to the infinite-order differential equation jSJ, let us now formally 
consider a set of truncated approximations of the BK equation: 

^Y^f^xp,',A'^L)^f-^P (4) 

where each truncated kernel XPna given by the expansion of the original kernel x around 70 (0 < 70 < 1) up to order 
P>2 : 

XPno i-dL) = E i-dL - 70)^ = Y.{-\YA,dl . (5) 

P and 7o are then two parameters which define a given approximation of the original kernel. Each particular kernel 
XP.70 is polynomial in dh and its coefficients can be fully computed from the original kernel x ■ they are given by 

^ ■ Y(*+PH7n) ■ 
-4p=E(-l) ■, I ^0- (6) 

i=0 ^ 

For the sake of simplicity, we do not write the subscript P, 70 of the coefficients Ap. We want to consider kernels 
XP,'yo that are somewhat close to x- To be more specific, we shall limit our study to kernels that are positive definite 
and feature Aq > 0, ^1 < and A2 > 0. We also require that, just as the original BK equation ((SJ, the equations 
Q admit asymptotic traveling-wave solutions. In other words, for a kernel XPno to be considered, the equation 
XP-iail) /1 — x'p -ygil) >0 should have a unique solution 7^''''°, which we shall call the critical anomalous dimension. 
We shall also denote the corresponding critical velocity v^''^" . 



For the kernels Xp,7cj fo'" which the expansion has been done around 7c, one has of course =7c and —Vc- 
Then to approximate the solutions of the BK equation (O by solutions of (^J , one considers kernels expanded around 
7o~7j, to insure that the speed of the resulting traveling wave is close to Vc- We also expect that a relevant truncation 
is obtained with a small value of P. 

Inserting JSJ into l@J, the truncated BK equation reads: 

p 

AqU-W"- dyU - Ai OlN + ^(-l)PAp dlM = . (7) 

p=2 

The advantage of this equation w.r.t. the original BK equation is that it can be exactly solved ^| in terms of 
traveling waves. Indeed, consider the anzatz 



with the scaling variable s given by 



M{L,Y) = Ao U{s) , (8) 



^^L-U + -^A,]Y (9) 



where A= ^'AqJA^ and c is a parameter related to the speed of the wave. The equation for the traveling wave U then 
becomes the following ordinary differential equation 

C/(l ~U) + U' + ^U" + C/(f) - . (10) 



, CP Aa 



Let us define 



v{c) = Ai+c^A^A2 (11) 



such that s~X{L — v{c)Y)/c and that u(c) is the speed of the traveling wave. For a given kernel XP,-yai the values of 
O — All ^ AqA2 for which equation ((Tmi has a solution define the possible values of the speed v{c) of the travehng-wave 
solution 

It is crucial to notice that, in this approach, c is a free and thus adjustable parameter. It means that in general, 
there is a continuum of solutions with different speed to the equation Q or equivalently 0. For instance, if one 
wants to describe asymptotic solutions, we are led to choose the value of c such that 

^;(c)=«f'-"'^Xp,^„(7f'^") if7^>7c 



(12) 

v{c) =XP.,ia{lt)llt ifO<7i<7c ■ 

Indeed, once the initial condition is fixed, the asymptotic speed is v^'^° or XP,7o(7i)/7i depending on the initial 
condition, as confirmed by numerical simulations. We shall use the freedom on c to also examine non-asymptotic 
properties. 

Equation (|l()|l provides a 1/cf expansion in with the 1/c term missing. The key point of the method is that 1/c is 
indeed a small parameter llH allowing to find an iterative solution 



h{s) = hois) + l/i2(.s) + J2 ^^(^) ^ J - ^(^) ■ (13) 

p>3 

Inserting it into (|10|) translates into a hierarchy of equations 



h'o + hl-l/A = 

h'.2 + 2hoh2 + h'^ = 

h'^ + 2hoh3- X^Ash'^'/Ao = 

h'i + 2hoh4 + hl + h'^ + X'^A4h'^"/Ao = 



(14) 



where we have written the equations up to 0{l/c^). An iterative solution can easily be found with initial conditions 
being appropriately chosen, i.e. ft.o(±oo) = ±i and hi^o{zLoo) = hi{0) = 0. Note that if h{s) is solution, also h{s + SQ) 
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FIG. 1: The function U{s) solution of the equation U{l-U)+U' + U" /c^ ^0 for c= 5/^6 with f/(-oo) = l and [/(0) = l/2. FuU 
line: exact solution HSU . Dashed line: solution 11611 containing only the first order. Dotted line (hardly distinguishable from 
the full line): solution p6^ containing only the first and second order. 

is solution; this freedom corresponds to the possibihty of arbitrarily redefining fco. The system is iteratively but fuUy 
solvable. One first solves the only non- linear equation of the hierarchy, obtaining /iq = ^ tanh(|). Using the property 



-rhnis) + 2hohn{s) 
as 



1 



cosh^(s/2) ds 



[cosh2(s/2) Kis) 



(15) 



all the other linear equations reduce to simple integrations of functions defined recursively from the hierarchy. The 
solution of the first three terms of the expansion is 



Uis) 



1 



1 



1 + e^ 



d + e'^ 



log 



4e* 



A3 A, 



c3 Ao {l + e'Y 



3?^ T + s 



(16) 



The first two terms of the expansion l|16(l (order l/c" and 1/c^) are universal: they do not depend neither on P nor 
on 7o. Indeed equation (@J) with any kernel XP,ia admits solutions whose first two terms (in the 1/c expansion) are 
those of Hlfci|) . The solutions for different kernels only differ through the definition ^ of the scaling variable s which 
depends on P and 70. In this sense, we obtain a universal parametric solution. 

Let us concentrate on the case P = 2. The equation (@J) reduces to the known FKPP equation and one obtains the 
corresponding traveling-wave equation: 



U{1 -U) + U' + —U" = . 



(17) 



One has ^2,70 {li)ht^Aohi+Ai+Mli and v^-'^" = Ai+2yy A0A2 and therefore, using (O one should take c = A/7i47i/A 
in the "pushed front" case and c = 2 in the "pulled front" case to reproduce the correct critical speed. 

In order to check the efficiency of the iterative solution, we use the existence of an analytic solution to equation 
(|17|l . for the particular value of c = 5/V6~2.04 : 



U{s) 



(V2 - 1) exp ( 



(18) 



This exact solution allows one to see that limiting the expansion after the second order is a very good approximation: 
in Fig. 1 we have plotted the exact solution p8|l along with the solution given by the expansion 1)16(1 limited after first 
and second order with c = 5/V6. It is very hard to distinguish the exact solution from the solution at second order. 



B. Comparison with numerical solutions of the BK equation 



In terms of physical variables, the solution (|16|l of equation ^ reads: 



Af{k,Y) 



Ao 



Ao 



fc2 



X/c 



A/c 



f^2 



x/c 



log 



X/c 



fc2 



UY) 



X/c 



(19) 



where Ql{Y) = kQ exp [ati(c)F] plays the role of the saturation scale and where we have restored the a dependence of 
Y. In the following we fix a = 0.2. The Ap coefficients are given by A=\/^o/^2 and c is chosen so that the traveling 
speed v{c) — Ai+c^/AqA2 takes the value of the critical speed v^'"^" in the "pulled front" case and XP,ioili) lli 
"pushed front" case, see H12|l . In the following, we shall only use the first two terms given in (|19|l which provide a 
good approximation of the full traveling-wave solution. 

Let us compare numerical solutions of the BK equation with the solution (|19|l for P — 2. For the numerical 
simulations, the initial condition we chose^ is Af{L,Yo) = e~'^'^ /{2^i) if L>0 and J^{L,Yo) — {l~jiL)/{2^i) if L<0. 
As already mentioned, only the large— L behavior of the initial condition is important to determine whether one lies 
in the pushed or pulled- front case. For the analytical solution (|19|l . the value of c that gives the critical speed is 
c — X/'^i+'ji/X in the pushed-front case and c = 2 in the pulled-front case for which we recall the value 7c = 0.6275. 

The comparison is shown in Fig. 2: the "pulled front" case is represented on the left plot for ^i — l and with the 
kernel expanded around 70 = 7c, the critical value. The "pushed front" case is represented on the right plot for 
7o = 7i = 0.5. One sees that as soon as the asymptotic traveling- wave regime is reached, there is a good description of 
the "interior" of the wave, as defined in Ref.'^. In both cases, one sees that equation (|19|l does not describe neither 
the saturated region (small k) nor the "leading-edge" (large k). Indeed, the exact scaling properties of traveling waves 
are mathematically expected only in the interior region 7|. 
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FIG. 2: The traveling wave N{L, Y) (fixed coupling) at large y as a function of L. The difi'erent values of Y are 0, 20, 40, 100. 
The full lines are numerical solutions of the BK equation and the dashed line are the scaling solutions 11911 for P = 2. Left plot, 
"pulled front" case (70 =7c, 7i = 1): the asymptotic traveling-wave regime is reached after some evolution with the critical 
speed «c = 4.883. Right plot, "pushed front" case (70 = 7i = 0.5 < 7c): the asymptotic regime is driven by the initial condition 
and the critical speed is X2,i/2(i/2)/(i/2) = 5.55. 



For our purpose, it is useful to investigate the non-asymptotic regime since it is of interest for phenonienological 
applications. In Fig. 3, we compare the numerical solutions of the BK equation with the solution p9|l . in the pulled 
front case, where we adjusted the parameter c to describe the speed of the wave v at moderate values of the rapidity 
Y up to 20. Again the agreement in the interior region is satisfactory and it shows that the solution we obtained can 
be used also in the physical range of non-asymptotic rapidities. 



^ This form ensures correct infrared and ultraviolet behaviors with an appropriate matching. The details of this matching are not relevant 
for the asymptotic properties. 
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FIG. 3: The traveling wave N{L, Y) (fixed coupling) at moderate y as a function of L. The different values of Y are 0, 2, 4, 20. 
The full lines are numerical solutions of the BK equation and the dashed line are the scaling solutions 1191 for P = 2, 70= 7c, 
and 7i = 1. The speed of the wave has been adjusted to the value v(c) —4.2. 



Let us finally point out the role of the kernel truncation. First, the truncation up to a finite P is mathematically 
necessary to apply our method. Indeed, when considering the full analytic expansion to define the coefficients Ai 
(i.e. P — > 00), they are infinite due to the 7 = singularity of the BFKL kernel. This explains why the solutions 
(|19|1 describe the interior region of the BK front but not the whole k range. For instance, one does not describe the 
saturated region formula (|19|l describes fronts which saturate at Aq while the BK front behaves as log 1/k when 
k is very small. A second remark is that we observed, e.g. by comparison of the results of the P = 2 and P = 4 
truncations, that a good description of the interior region of the wave is already obtained for P = 2, i.e. the diffusive 
approximation of the BFKL kernel. In this paper, we shall thus stick to P = 2. 



III. THE RUNNING-COUPLING CASE 



Let us now perform the analysis in the running-coupling case with 

Nr^ ^.2A 1 , lliV,-2iV/ 

— ' - 1 - 

'■QCD 



and Aqcd = 200 MeV. We thus write the following form of the BK equation with running coupling constant and 
leading-order BFKL kernel ((SJ: 

bLdyM^xi-dL)^-^^ ■ (21) 



It has been proved |2| that this equation admits asymptotic traveling-wave solutions of the form N [L — ^Ji^IvlVyY) 
where as before v — Vc in the "pulled front" case and v — xili)/ji in the "pushed front" case. 

Let us consider the same truncated kernels than in the previous section. That transforms the equation into 

p 

AoN-M^-bL dyM - Ai DlM + Y^{-lYAp dlN = (22) 
and we consider again the following form 

U{L,Y)=A^U{s). (23) 

Our strategy is to reduce the problem to the one we solved in the previous section through an appropriate change of 
variables. We require that the universal terms in the corresponding equation for U are the same as before: C/(l— [/)+[/'. 



Postulating s = L (piY/L"^) determines a solution for and leads to the scaling variable 



where c is a free parameter. Putting this into \22\ gives 

'■I , 

All Ao 



[7(1 _[/)+[/'+ ^ ( 4^ t/(p) + o 

p=2 



(24) 



(25) 



The leading-order terms in 1/c are of the same generic form as equation (jlOjl and thus one can apply the previous 
method. In this running-coupling case, the 0(1/5) terms that we shall neglect contain scaling-violation corrections 
decreasing as 1/L. Using the method described in the previous Section, one obtains the following solution 



Ui~s) 



AoY A2 



log 



4e^ 



(26) 



where the expansion parameter is Aq/Ai. As before, we shall only consider two first terms of the expansion. In order 
to describe solution to the BK equation (|21|l , one has now to select the speed of the wave by fixing the parameter c. 
Writing at large Y 



Ao 
Ai 



Aoc 



V-'^AiY 



and matching with the asymptotic form JV{L— y/ {2v/b)Y), one gets 

bAl 



A'o ■ 



(27) 



(28) 



Adjusting c to obtain a given speed, we check that 1/c is indeed small. Defining Ql{Y) — Aq^j exp{^/(2Mc)/l)yY) 
that plays the role of the saturation scale, the scaling variable reads in terms of physical variables 



QCD 



\ 



blog^ 



fc2 



^^QCD 



A\ 



■log' 



Formula along with 



Mik.Y) 



Ao 



KM 
A\ 



(1- 



log 



4e^ 



(29) 



(30) 



is the traveling-wave solution of the truncated BK equation with running coupling H21(l at C(l/c). 

Let us compare numerical solutions of the full BK equation l|21|) with that parametrization. For the numerical 
simulations, the initial condition is the same as in the fixed-coupling case. The problem of the Landau pole of the 
running coupling is dealt with by introducing a regulator: log(A:2/AQ(^^) ^ log(K-|-/c2/AQ(^^) with k = 3. As well as 
the details of the initial condition, this does not modify the traveling-wave pattern. 

The comparison is represented in Fig. 4 where the "pulled front" case is considered. The values of the coefficients 
Ao = 9.55, Ai = —25.6 and A2 = 24.3 have been fixed by considering P — 2 and 70 = 7c. We use the free parameter 
c to adjust the speed of the wave. On the left plot, high values of the rapidity Y — 0, 10, 20, 100 are considered 
while on the right plot, moderate values F = 0, 2, 4, 20 are represented. In both cases, in order to describe the 
numerical solutions, one has to adjust the actual speed of the wave v{c) to the value of 3.1 for the left plot and 2.3 for 
the right plot which are both significantly smaller than the critical value Dc — 4.88. It confirms that the establishment 
of the asymptotic regime is much slower in the running-coupling case than in the fixed coupling case Indeed, 
the theoretical sub-asymptotic corrections to the critical speed have been computed and predict such a behavior. 
There is a good description of the interior of the wave for the different ranges of rapidity, provided we adjust the 
speed of the wave. As expected, the parametrization does not describe the saturated and the leading-edge regions. 
Note that the scaling region at low values of Y is less extended in the running coupling case (Fig. 4b) than for fixed 
coupling (Fig. 3), this could be due to the scaling-violation terms that we had to neglect in 1)25(1 . However it is clear 
that an approximative traveling-wave pattern already emerges at moderate rapidities. 



-10 10 20 30 40 -10 -5 5 10 15 20 

log(*:2) log(k^) 

FIG. 4: The traveling wave N{L, Y) with running coupling as a function of L. The full lines are numerical solutions of the 
BK equation with running-coupling and the dashed lines are the parametrization l|29|l and Left plot: high-energy regime 

y = 0, 10, 20, 100 and ^(c) =3.1. Right plot: moderate energy regime Y = 0, 2, 4, 20 and w(c) = 2.3. The initial condition is 
chosen in the pulled front regime. 



IV. APPLICATION TO DEEP INELASTIC SCATTERING 



We shall now draw a link between traveling waves and phenomenology. The observation of geometric scaling asks 
the question whether it can be viewed as a consequence of QCD traveling waves. The known problem is that the 
range of rapidity for which asymptotic predictions of QCD non-linear evolution equations are available remains far 
from the moderate rapidity range accessible to experiments. Hence the phenomenological interest of our method is 
to open a way of confronting theory and experiment in some physical kinematic range. 

The formula relating the dipole scattering amplitude in coordinate space N{r, b, Y) to the unintegrated gluon 
distribution of the target f{Y, fc^) is (see e.g. jl^): 

J dH Nir, h,Y)^^Jf f{Y, - Jo(fc|r|)) . (31) 

In the b— independent approximation, one writes / cPb N{r, b, Y) = ■KR^N{r, Y) where Rp is the radius of the proton 
target. Through formulae and H31|l a link is established between M{k, Y) and /(Y, k"^). One gets 

Finally, we use the relation 

f{Y,k') = ^ xgix,k') , (33) 
where F = log(l/a;) between the unintegrated gluon distribution f(Y, k"^) and the gluon distribution xg(x, k"^) to obtain 

Note that formula H33|l is a well-known approximation which may only be valid at small values of x. It is useful here 
as it makes the analysis simpler. For a more detailed study which is beyond the scope of this work, one could consider 
more advanced prescriptions |14| . 

In order to compare (|34|l to our predictions, we use the running-coupling case which is known to lead to traveling 
waves with a speed compatible with phenomenology |l5l |. 

We performed a fit of the parametrization (|29|l and (|30l) on the amplitude (|34|l obtained from the Martin-Roberts- 
Stirling-Thorne (MRST) gluon distribution 16] . The domain we considered is x < 10~^ and ^l^U{k, Y) > 0.01 which 
is the region where traveling-wave patterns can be tested. Note that we considered only the first term of (|30l) which 
depends on the parameters Ag, Ai, and v the speed of the wave. We discarded the second term of since we found 
that the parameter A2 stays undetermined by the fit in the kinematical region we considered. In Fig. 5, we show the 
result of the fit, the obtained parameters being Aq — IT.I^ Ai = — 15.8, and v = 1.7Q. 

These results call for comments. 
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FIG. 5: The dipole amplitude N{k, Y) as a function of L for different values of y = 4, 5, 12. The full lines are obtained from 
the MRST gluon distribution via 134^ and the dashed lines represent the parametrization H29|l and lj,SO|l . Only one term is kept 
in CT . 

i) The dipole amplitude obtained from the MRST distribution (full lines on Fig. 5) displays a structure compat- 
ible with an evolution towards traveling waves, namely a steeper slope at small rapidities which evolves toward a 
less steep and regular pattern at higher rapidities. 

ii) Fig. 5 shows a compatibility with the parametric form f (|29|l . (|30|l 'l at high values of Y and moderate values of k^. 

iii) On the theoretical ground, it is worth comparing the fitted parameters to various predictions. For the truncated 
BK kernel with P — 2 and op = 7c, one finds Ao — 10, Ai = —25, and the speed 2.3 < v <3 when one includes 
pre-asymptotic corrections 0. With P = 4 and 70 = 1/2, one finds Aq — 11, Ai = —49. The obtained parameters 
are qualitatively closer to the diffusive approximation (P = 2) of the BK equation but they are quantitatively 
different. This seems to suggest the presence of next-leading or higher-order corrections to the QCD kernel. This 
deserves further study. 

iv) Note that we obtain l^oMil ~ 1 for the observed parameters which could invalidate the expansion H26|l . We 
checked both numerically and from the analytical form that the second and higher order terms stay sufficiently 
small to be still neglected. 

V. CONCLUSION 

Let us summarize the main results of our study. 

- We found iterative traveling-wave solutions to non-linear BK evolution equations obtained by finite truncation of 
the BFKL kernel for both fixed and running coupling constant. 

- These solutions exhibit universality properties, the first two dominant terms of the iteration have a parametric form 
independent of the truncation of the kernel and of the fixed or running coupling cases. 

- The scaling variable is found to be a combination of L = log{k^) and rapidity Y which is given by formula © (resp. 
(|24|l ') for fixed (resp. running) coupling. 

- The obtained traveling-wave solutions match with the asymptotic solutions of the BK equation. This was verified 
by analytical and numerical checks. The remarkable new property is that they also match with the non-asymptotic 
behavior of the interior of the wave, provided an adjustment of the speed of the wave which is a free parameter in the 
iterative approach. 

- As an application of the method and its validity at non-asymptotic energies, we considered the dipole amplitude in 
momentum space Af{k, Y) obtained from the MRST parametrization of the gluon distribution. We found evidence 
for an evolution pattern compatible with the formation of traveling waves. The obtained dipole amplitude is well 
described by the universal parametric form (|30|l . The obtained parameters seem to point towards traveling- wave 
solutions for a BFKL kernel modified by higher orders and/or non-perturbative corrections. 

Our results help establishing a stronger link between the theory of saturation and the phenomenological observation 



of geometric scaling Moreover, simple analytical parametrizations of the solutions of the BK equation in momen- 
tum space could help extending the present phenomenology |17| to a broader range of kinematics or observables. It 
would be interesting to investigate where one could distinguish between a DGLAP structure and traveling waves. On 
a more theoretical level, it seems feasible to extend the method to QCD evolution equations beyond the "mean field" 
BK equation In particular the extension to the stochastic versions of QCD evolution equations |0| would help 
our understanding of QCD in the high-energy limit. 
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